o 

o 

(N 

Ph. 
< 






o 
m 



RBRC-742 

On-shell A/ = 3/2 kaon weak matrix elements with non-zero total 

momentum 

Takeshi Yamazaki *^'^ 

(RBC and UKQCD Collaborations) 

^Physics Department, University of Connecticut, 

Storrs, Connecticut 06269-3046, USA 

^RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973 

(Dated: April 14, 2009) 

Abstract 

We present our results for the on-shell A/ = 3/2 kaon decay matrix elements using domain wall 
fermions and the DBW2 gauge action at one coarse lattice spacing corresponding to a~^ = 1.31 GeV 
in the quenched approximation. The on-shell matrix elements are evaluated in two different frames: 
the center-of-mass frame and non-zero total-momentum frame. We employ the formula proposed by 



Cn ' Lellouch and Liischer in the center-of-mass frame, and its extension for non-zero total momentum 

> 



frame to extract the infinite volume, on-shell, center-of-mass frame decay amplitudes. We determine 
the decay amplitude at the physical pion mass and momentum from the chiral extrapolation and 



C"^ I an interpolation of the relative momentum using the results calculated in the two frames. We have 

g ■ obtained Re^a = 1.66(23)(to3)(lf ) x 10"^ GeV and Im^s = -1.181(26)(tJ^^)(1^4) x lO-^^ GeV 

at the physical point, using the data at the relatively large pion mass, m-„: > 0.35 GeV. The first 
r> I error is statistic, and the second and third are systematic. The second error is estimated with 

several fits of the chiral extrapolation including (quenched) chiral perturbation formula at next to 
leading order using only lighter pion masses. The third one is estimated with an analysis using the 
lattice dispersion relation. The result of Re^2 is reasonably consistent with experiment. 

PACS numbers: 11.15.Ha, 12.38.Aw, 12.38.-t 12.38. Gc 



present address: Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, Kyoto 606-8502, 
Japan 



I. INTRODUCTION 



There are longstanding problems in non-leptonic kaon decays, such as, A/ = 1/2 selection 
rule and the ratio of the direct and indirect CP violation parameters e' je [l|, l^] , which need 
to be resolved for precision tests of the standard model. The major difficulty is to evaluate 
the non-perturbative strong interaction effect. Although such non-perturbative effect in 
principle can be quantified by lattice QCD, there are many technical difficulties to calculate 
the if — > vrvr decay process directly on the lattice. This is because there are problems to 
deal with the two-pion state in finite volume [3|. To avoid the difficulties, Bernard et al. 
proposed one possible approach, called as the indirect method, [^ where K ^> ixix process 
is reduced to if ^ vr and if ^ processes through chiral perturbation theory (ChPT). So 
far several groups have reported results obtained by this method p, IgL LZI, Isl, la, [lO|, lll |. 

Especially, it is worth mentioning that CP-PACS [lOi] and RBC [ll| collaborations have 
independently calculated the full weak matrix elements of the if ^ tttt decay with the 



indirect method by using domain wall fermion action p^, Il3|, IM], which has good chiral 
symmetry on the lattice. Their final results of e' je^ however, have the opposite sign to that 
of the experiment. Indeed, in their calculations there are many systematic uncertainties since 
those calculations have been performed at a single finite lattice spacing using the indirect 
method based on tree level ChPT within the quenched approximation. The indirect method 
might cause the larger systematic error than other sources, i.e., the quenching effect and the 
scaling violation, because the final state interaction of the two-pion is expected to play an 
important role in this decay process. Therefore we have to directly calculate the scattering 
effect of the two-pion state on the lattice to eliminate this particular systematic error. For 
this purpose we attempt to carry out a calculation with direct method, where the two-pion 
state is properly treated on the lattice, in the Ai = 3/2 K —* tctt decay process. 

There are two main difficulties for the direct method. The first problem is that it is 
hard to extract two-pion state with relatively large momentum, e.g., about 200 MeV, on 
the lattice by a traditional single exponential analysis. The problem was pointed out by 
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IGj calculated the two-pion 



Maiani and Testa [3|. To avoid this problem two groups 
at rest, but allowing a non-zero energy transfer in the weak operator. We however cannot 
obtain physical amplitudes from the calculation, unless we use some effective theory, such 
as ChPT, to extrapolate unphysical amplitudes to physical one. 



There are several ideas for solving the problem. One of the ideas is to employ a proper 
projection of the K ^ -n-n four-point functions [l7|. In this approach we need compli- 
cated calculations and analyses, e.g., diagonalization of a matrix of the two-pion correlation 
functions 18|, to treat the two-pion state with non-zero relative momentum on the lattice. 



A simpler idea, where such complicated analyses are not required, is to prohibit the 
two-pion state with zero momentum. Recently Kim |19l . |20(] reported an exploratory study 
with H-parity (anti-periodic) boundary conditions in the spatial direction. The method 
works to extract the two-pion state with non-zero relative momentum from the ground state 
contribution of correlation functions, because the two-pion state with zero momentum is 
prohibited by the boundary condition. 

Alternatively, if we perform the calculation in non-zero total momentum (Lab) frame, 
we will be able to forbid the zero momentum two-pion state, which certainly appears in 
center-of-mass (CM) frame. In the simplest Lab frame the ground state of the two-pion is 
|7r(0)7r(P)) with P being non-zero total momentum, which is related to the two-pion state 
with the non-zero relative momentum in the CM frame. Therefore, we can extract the two- 
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23| , as well 



pion state with non-zero momentum from the ground state contributions |2lL 

as in the H-parity boundary case. In this work we mainly employ the latter, namely the 

Lab frame method. 

The other difficulty of the direct method is the finite volume correction due to the presence 
of the two-pion interaction. We have to estimate this finite volume effect to obtain desired 
matrix elements in the infinite volume, because the effect is much larger than that of one- 
particle state. Lellouch and Liischer (LL) [2J] suggested a solution of this difficulty, which is a 
relation of the on-shell, CM frame decay amplitudes in finite and infinite volumes. However, 
their derived relation is valid only in the CM frame with periodic boundary condition in the 
spatial direction, so that we need a modified formula when we utilize H-parity boundary 
condition 20| or Lab frame method. Recently, two groups, Kim et al. 25|| and Christ 
et al. 261], generalized the formula that is an extension of the LL formula for the Lab 
frame calculation. The generalized formula is based on the finite volume method to extract 
the scattering phase shift from the Lab frame calculation 23] derived by Rummukainen 
and Gottlieb. Here we attempt to apply this generalized formula to the calculation of the 
A/ = 3/2 kaon weak matrix elements with domain wall fermions and quenched DBW2 



gauge action [271, l28| at a single coarse lattice spacing. 



It is worth noting that there are some (dis) advantages of the two methods, the Lab frame 
and H-parity boundary condition calculations. In the H-parity method, the statistical error 
of the matrix elements is smaller than those in the Lab calculation when same simulation 
parameters, e.g., lattice size and quark mass, are utilized. This is because the unit of the 
momentum in the H-parity method is a half of that in the Lab frame calculation due to the 
anti-periodic boundary condition in the spatial direction. However, the H-parity method is 
not effective in the calculation of A/ = 1/2 K ^ tttt decay, because under the H-p arity 



boundary condition only vr^ states satisfy the anti-periodic boundary condition |19l . |29 |. 
The vr'^ state does not satisfy the boundary condition, unless we introduce another discrete 



boundary condition, such as G-parity boundary condition J30| . On the other hand, the Lab 
frame method does not break the isospin symmetry, so that we can apply the method to the 
calculation of AJ = 1/2 channel as well as that of AJ = 3/2. 

The organization of the article is as follows. In Sec. [Ill we give a brief explanation of 
the generalized formula as well as the original LL formula. We also explain the calculation 
method of the correlation functions, from which we evaluate the decay amplitudes. The 
parameters of our simulation are given in Sec. IIIII We show the result of the 1 = 2 two-pion 
scattering length and phase shift in Sec. IIV[ Then we present the results for the off-shell 
and on-shell decay amplitudes, weak matrix elements, R.eA2, and Imy42. Finally, we br iefly 



summarize this work in Sec. |Vl Preliminary result of this work was presented in Refs. 3ll.l32| 



II. METHODS 



A. LL formula 



Let us briefly review the method suggested by Lellouch and Liischer [24] (LL). They 
derived a formula, which connects a CM decay amplitude \A\ deflned in infinite volume to 
the one given on finite volume \M\ as, 

where Ej^^^ is the two-pion energy in the CM frame and the scattering phase shift 6 is 
responsible for the tttt final state interactions. The relative momentum of the two pions p is 
determined from measured i?^^, 

p' = El/A - ml (2) 



and its normalized momentum q is defined by 

q' = ipL/27rY. (3) 



with the spatial extent L. The function 0(g), derived by Liischer [33|, is defined by 



g7r3/2 
tan0(g) = - (4) 



V47r f-', n^ - g^ 



where 



Using the function 0(g), we can also determine the scattering phase shift at the momentum 
p through the relation 

S{p) =mr-(j){q), (6) 

where n is an integer, and 0(0) = 0. Note that the formula eq. ([1]) is valid only for on-shell 
decay amplitude, i.e., E^^t, = uik- 

B. Extended LL formula for non-zero total momentum frame 



Next, we describe the extended LL formula for the the non-zero tota! 



momentum frame 



26| have derived a 



(denoted in the following by the Lab frame). Recently two groups [25, 

formula, which connects a CM decay amplitude |y4| in the infinite volume to a Lab frame 

decay amplitude \M^\ given on finite volume 

\A\^ = 8.7^ {^y{p^ + ^^} iM-r, (7) 

where 7 is a boost factor given by 

7 = E^JE^^, (8) 

with E^^ being the Lab frame value, which is determined by two-pion energy with the total 
momentum P. The total CM frame energy of two-pion states -Etttt is evaluated through the 
energy-momentum conservation, 

EL = (E^r - P'. (9) 

In eq.([7]), p and q are defined as the same in eqs.([2]) and ([3]), except the determination 
of E'^TT, which is given by eq.([9]). In the following the momentum denoted by the capital 



letter like P represents the total monientum, while the small one p represents the relative 
momentum, unless explicitly indicated otherwise. In addition quantities with a superscript 
P, e.g. E^^, indicate the Lab frame values. The function 4>p{q), derived by Rummukainen 



and Gottlieb 22j, is given by 



where 



tan0p(g) = - ^7^, (10) 

^oo(i;r;7) 
^°°^^' ^'' ^^^V^ ^, n? + ni + 7-^(^3 + 1/2)2 -g^' (11) 

in the case for P = (0,0,27r/L). The formula given in eq. ([7]) is valid only for the case 
E^^ = E^ where E^ is kaon energy with the momentum P, in other words, E^^^^ = thk 
as in the LL formula given in eq.([l]). When we set P^ = and 7 = 1, the formula eq.(I7]) 
reproduces the original LL formula eq. ([I]). Of course, we can determine the scattering phase 
shift at the momentum p by using (ppiq) as 

S{p)=n7r-<Pp{q), (12) 

where n is an integer, and (f)p{0) = as the same in eq.©. 

C. Calculation of K ^ tttt four-point function 

We calculate four-point functions for the A/ = 3/2 K -^ tttt decay in zero (CM) and 
non-zero (Lab) total momentum frames. In actual simulations, we employ P = (0, 0, 27c/L) 
for the non-zero total momentum. The four-point functions in the CM and Lab frames are 
defined by 

Gi{t,U,tK) = (0|7r+7r°(0,tOOf/'(t)[ir+(0,tx)]^|0), (13) 

G^it,U,tK) = (0|7r+7r°(P,t^)of (t)[K+(P,t^)]t|o), (14) 

where the two-pion operator in Gf is averaged with the following two types of products of 
pion operators in order to extract the 1 = 2 part of the two-pion state, 

7r+7rO(P, U) = (7r+(P, t^)n%6, U) + 7r+(0, t^)7r%P, U))/2. (15) 



Oj are lattice operators entering AJ = 3/2 weak decays: 



of = ^|(s"(f)rf"(f))L[(M^(f)M'(f))L-(/(x)rf^(f))J + {r{x)u''{x))L(u\x)d\x)) 



(16) 



0%' = 5^{(s"(f)rf"(x))L[(M'(f)M^(f))ij-(/(f)rf'(f))J +irix)u^ix))Liv'ix)d\x))R 



(17) 



,3/2 



0;;'88 = $^](s''(f)f/'(f))L (tl'(f)w'^(f))iJ-(/(f)t/'^(f))K +{r{x)u\x))L{v'{x)d^{x))j 



where {qq)L = ^7At(l ~ 75)^') (j1'1)r = ^7/^(1 + 75)'?; and a, 5 are color indices. The kaon weak 
decay operators are classified into the (8,1), (27,1), and (8,8) representations of SU{3)l ® 
SU{3)r, but the A/ = 3/2 part has only the (27,1) and (8,8) representations. O27 and Ogg 
are the operators in the (27,1) and (8,8) representations with / = 3/2, respectively. O^gg 
equals Ogg with its color summation changed to cross the two currents. 

We employ the momentum projection source for the quark operator with the Coulomb 
gauge fixing to obtain better overlap to a state with each momentum. The pion operators 
with the source are constructed as, 



TT^ilt) = [J2dix,t)e-''-'Mj2<y,t)e 



-ik2-y] 



(19) 



^\k,t) = [^^(f,t)e-^'-%5E«(y,t)e-^'-^1-[X^rf(f,t)e-^'=-^l75E^(^>^)e"^'"'] I /v'2, 

\ X y X y / 

(20) 

with k = ki + k2 = P OT 0, where the momentum k and ki represent the momentum of the 
pion and each quark, respectively. In the zero momentum case, these operators are nothing 
but the wall source operator. The K~^ operator with each momentum is calculated in the 
same way with changing c? to s in eq. (ll9l) . 

The four-point function in the Lab frame is much noisier than that in the CM frame. In 
order to improve the statistics, we calculate the four-point function Gf (t, t,r,^x) with two 
possible momentum insertions, ki = 0,k2 = P and ki = P,k2 = in eqs. OlQI) and (l20l) for 
the pion operators, and then we average them on each configuration. On the other hand, 
we fix the momentum of the kaon operator in the four-point function as ki = P, /c2 = 0. 



We also calculate four-point function for the 1 = 2 two-pion and the two-point function 
for the kaon and pion with zero momentum, 

GUt,U) = (0|7r+7r+(0,t)[7r+7r+(0,t^)]t|0), (21) 

GK{t,tK) = {0\K+{6,t)[K+{6,tKW\0), (22) 

G^{t,t^) = (0|7r+(0,t)[7r+(0,t^)]t|0), (23) 

and with non-zero total momentum, 

G^^{t,Q = (0|7r+7r+(P,t)[7r+7r+(P,t.)]t|0), (24) 

G^(t,tx) = {0\K+{P,t)[K^{P,tK)]^0), (25) 

G^{t,U) = (0|vr+(P,t)[7r+(P,t^)]t|0). (26) 

We employ the operator with ki = 0,k2 = P in eq.( IT9l) for G^it^tT,) and G^^itjt-,,), while 

with ki = P,k2 = in the kaon operator of G^{t,tK)- At the sink t of Gt„t ,Gj^ , and 

(p) 
Gn , we use the same operators as the source operator, and also the momentum projected 

operator for the meson field. The 7i~^{k,t) field with the latter operator is given by 

5^rf(f,t)75n(f,t)e-*^^ (27) 

X 

with fc = P or 0. In the zero momentum case, the operator corresponds to the point sink 
operator of the pion. We shall call the symmetric correlator under exchange of the source 
and sink operators as "wall" sink correlator, while the one calculated by the different source 
and sink operators as "point" sink correlator. We determine the energy for each ground 
state from the point sink correlator, and the amplitude from the wall sink operator. 

We use linear combinations of quark propagators with the periodic and anti-periodic 
boundary conditions in the temporal direction to have a quark propagator with 2T period- 
icity in the time direction. 

III. SIMULATION PARAMETERS 

Our simulation is carried out in quenched lattice QCD employing a renormalization group 
improved gauge action for gluons, 

f3 



Sg[U] 3 



:i-8ci) J2 p[uu,.+ci Yl RPUt^^ 



8 



(28) 



where P[U]x,iiu and R[U]x,iiu represent the real part of the trace of the path ordered product 
of hnks around the 1x1 plaquette and 1x2 rectangle, respectively, in the /i, z/ plane at 
the point x and (3 = Q/g"^ with g being the bare coupling constant. For the DBW2 gauge 



action 27|, |28|], the coefficient Ci is chosen to be —1.4069, using a renormalization group flow 



for lattices with a ~ 2 GeV 



27 



Q. 



34| 



Gauge conflguration was previously generated 
at /? = 0.87 with the heat bath algorithm and the over-relaxation algorithm mixed in the 
ratio 1:4. The combination is called a sweep and physical quantities are measured every 200 
sweeps. 



We employ the domain wall fermion action [l^, |l3|, \l^ with the domain wall height 
M = 1.8 and the flfth dimension length L^ = 12. The residual mass is reasonably small 
in these parameters, rrires = 0.00125(3) 35|. Our conventions for the domain wall fermion 



operator are given in Ref. [36|. The inverse lattice spacing is 1.31(4) GeV [34l] determined 
by the p meson mass. The lattice size is L^ x T = 16'^ x 32, where the physical spatial extent 
corresponds to about 2.4 fm. 

We flx the two-pion operator at t^ = 0, while we employ three source points tx = 
16, 20 and 25 for the kaon operator to investigate tx dependence of the statistical er- 
ror of the Lab frame decay amplitude, and to check the consistency of these results. 
We employ four u, d quark masses, m^ = 0.015, 0.03, 0.04 and 0.05 corresponding to 
m^ = 0.354(2), 0.477(2), 0.545(2) and 0.606(2) GeV, for the chiral extrapolation of the decay 
amplitudes. 

To utilize the LL formula eq.([I]) and the extended formula eq.(I7]), it is important to obtain 
a decay amplitude at on-shell, where the energies of the initial and flnal states are equal. In 
order to obtain the decay amplitude, we vary the kaon energy with several strange quark 
masses at a flxed light quark mass. Then we carry out an interpolation of the amplitudes 
to the on-shell point. For the interpolations, six strange quark masses are employed, rris = 
0.12, 0.18, 0.24, 0.28, 0.35 and 0.44, except the lightest m„ for the t^ = 16 and 20 cases where 
we use the three lighter rris- We will see that the three strange quark masses are enough for 
the interpolation of the case in a later section. 

In our simulation we employ relatively heavy strange quark masses to obtain the kaon 
energy closely satisfying to the on-shell condition, Ej^^^ ~ itik- This might cause a systematic 
error, because lower modes of the hermitian dirac operator would not be well separated to 
left- and right-hand in a heavy quark mass as discussed in Ref. [32]. The systematic error 



can be removed by simulations with lighter pion mass, larger volume, and finer lattice 



spacing 37|. However, we will not estimate such effect in this work, which may be accepted 
as an exploratory work. Thus, it should be reminded that our result might suffer from above 
mentioned systematic error at least at the heavier strange quark masses. 

The LL formula and its extension requires a relatively large spatial volume. This require- 
ment stems from both the on-shell condition and size of the scattering range R. R is defined 
by 

V{r) = at |r=l > i?, (29) 

where V{r) is the effective potential of the two-pion scattering at the relative coordinate 
r of the two pions. The former is much essential here. The on-shell condition requires a 
large spatial extent, e.g., L = 6(3.2) fm for the CM(Lab) frame, which is evaluated at the 
physical pion and kaon masses. The required spatial size, indeed, changes in accordance 
with simulated masses of the pion and kaon states, as L^ = ATr'^/{m'j^/A — m^) for the CM 
frame and LF' = 47r^/(?Ti|^/4 — ml) ■ m^/m\ for the Lab frame. Thus, it is possible to fulfill 
the on-shell condition on a smaller volume than the one required at the physical point. In 
this work we tune itik at the simulated m^ and the fixed L to satisfy the on-shell condition 
in each frame as we explained above. 

Another constraint on the spatial volume is determined from the scattering range R, 



because R < L/2 is an important assumption [22|, |33|] to derive the two relations eqs.([6]) 
and flT2l) . where we can evaluate the scattering phase through the two-pion energy on finite 
volume. The LL formula and its extension are based on those relations. If the assumption is 
not satisfied, the scattering effect is contaminated by unwanted finite spatial size effect, in 
other words, the effective potential is distorted by the boundary condition. The determina- 



tions of R were carried out numerically in both the CM [38| and Lab [39] frame calculations 
in the 1 = 2 two-pion channel. These references reported that the required spatial extent 
was estimated as L = 2.4-3.2 fm in the range of ?7i^ = 0.42-0.86 GeV. The spatial volume 
~ (2.4 fm)^ in this calculation is not fully justified in this sense. Thus, we simply assume 
that effects stemming from the distortion of the effective potential are small or negligible in 
our simulation. 

We use 252 configurations in the case of tx = 16 and 20, except at the lightest pion mass 
where 371 configurations are employed in order to improve statistics, while we employ 100 
configurations in the case of tx = 25 at all the quark masses. The numbers of the gauge 

10 



configuration used in our simulations are summarized in table [H 

In the following analysis we evaluate the statistical error of all the measured quantities 
by the single elimination jackknife method. For chiral extrapolations in the tx = 16 and 
20 cases, we employ the modified jackknife method [401] to take into account the different 
numbers of the configurations between in the lightest and other heavier u, d quark masses. 



IV. RESULTS 



A. Physical quantities for / = 2 vrvr scattering 



The calculation of the two-pion scattering is important not only to employ the LL type 
methods, but also to understand hadronic scattering from lattice QCD. So far many groups 



ca 
46 



culated the scattering length of the S-wave I = 2 tttt channel 



m, 



38 



41 



42 



43 



44 



with the finite volume methoc 



50. 


51 



33 



45|, 



49| . Some of the works obtained the result 



in the continuum limit |23l . 



53l |. Recently, the result using domain wall valence 
quark on the 2+1 flavor improved staggered sea quark was reported by Beane et al. 5J, |55 |. 
There are a few papers to calculate the I = 2 tttt scattering phase shift, with the finite 

m 



19 



20 



48] and 



volume method and its extension [22], at only a single lattice spacing 
in the continuum limit with quenched approximation 53| and two-fiavor dynamical quark 
effect 23|. It is worth remarking that there is the recent work for the I = 1 tttt scattering 
phase shift with p meson resonance 56| . In this section we present our results of the 1 = 2 
TTTT scattering length and scattering phase shift. 



1. Scattering length 

We evaluate the scattering phase shift 6{p) through the finite volume method eq.([6]) of 
Liischer |33| in the CM calculation. In the CM calculation we extract the two-pion energy 
£^,r7r to fit the point sink two-pion correlator using the fit form |20| . 



A.(« 



hz-^T^t 



+ e 



-E^^{2T-t) 



+ C) 



(30) 



We should notice that we employ linear combinations of quark propagators with the periodic 
and anti-periodic boundary conditions in the temporal direction, so that the periodicity of 
the correlators is 2T. The fitting analysis is carried out with three parameters A, C, and 
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i^TTTT in the fit range of t = 6-31. The constant C stems from two pions propagating opposite 
way in the temporal direction due to the periodic boundary condition. 

Using the finite volume method eq.(jn]), we obtain the scattering phase shift from the 
relative momentum which is determined from E^^^^ through eq.([2]). The results for -E^tt, p, 
and 5{p) in the CM calculation are summarized in table [ITl 

We estimate the scattering length oq with the measured 5{p) by using an assumption, 

a„ « M, (31) 

P 

where the scattering length is defined by 

ao = lim^. (32) 

This assumption is valid in our CM calculation, because "p is small enough. The measured 
value of aQ/rriT, is presented in table [TTl Figure [1] shows that the measured value of aQ/m.,,, 
denoted by "CM analysis", has an appreciable curvature for the pion mass squared. The 
scattering length ao/m.,^ at the physical pion mass, m.,^ = 140 MeV, is determined by a chiral 
extrapolation with the fit form, 

A + Bml + Cml, (33) 

whose result is summarized in table IIIII In Fig. [H the result at the physical pion mass is 



compared with the experiment |57i] and a prediction of ChPT 58 1 

We also fit the data using the prediction of the next to leading order (NLO) ChPT |59| . 

1 



lixP 



^-^.U^^)--M^ 



(34) 



where /i is scale, q = 7/2, and l^ji) is a low energy constant. / is the pion decay constant at 
the chiral limit. The scale is fixed at /i = 1 GeV for simplicity. We cannot obtain a reasonable 
X^/d.o.f. from a fit with one free parameter /(yu), where we use / = 0.133 GeV [35]. As 
summarized in table HVl a three-parameter fit with /, /(/i), and q, gives a reasonable value 
of x^/d.o.f. While in the fit / is consistent with the one at the chiral limit |35|, the coefficient 
Q differs from the prediction of ChPT. Although the quenched ChPT formula 60|, l6l| is also 
available, the quality of the fit is similar to what we obtained with eq. (IM|) . At the physical 
pion mass the ChPT fit result agrees with the simple polynomial fit in the above, so that 
we choose the polynomial one in the following analysis. 
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2. Scattering phase shift 

The two-pion energy in the Lab frame is noisier than the one obtained in the CM cal- 
culation, so that we determine the two-pion energy in the Lab frame E^^ in the following 
way to reduce the statistical error. First, we fit the point sink two-pion correlator with the 



non-zero total momentum using the fit form 



2l| 



where E^^ is the measured value of the single pion energy with the momentum. The terms 
multiplied by C come from two pions propagating in different temporal direction due to 
periodic boundary condition as in eq. fl30|) . The fit is carried out with three parameters A, 
C, and W in the range of t = 6-3L Using the parameter W, we determine the energy shift 
^E^^ = W — {E^ + m,r) from the non-interacting two-pion energy. Finally we reconstruct 
E^^ using m^r, AE^^, and the total momentum squared P^ 

E^^ = AEf, + m. + ./ml + P^. (36) 

We determine the CM frame, two-pion energy E^^t, from E^^ with eq. ([9]). We also 
determine the relative momentum p from E-j^^^ as in the CM case through eq.([2]), and then 
obtain the scattering phase shift by the extension of the finite volume method to the Lab 



frame, eq.(IT^ [22|. The results for E^^, ^E^-ni -^tttt, 7 = E^^/E^,^, p, and 5{p) in the Lab 
calculation are summarized in table |Vl 

To investigate the momentum dependence of the scattering phase shift, we use the results 



obtained from both the frames. We define "scattering amplitude" T{p) as in Refs. 23|, |48|, 



53 



62| 



Tip) = '-^^^M . ^, (37) 

which is normalized by aomT^ at the zero momentum. The results of T{p) in each frame are 
tabulated in tables [TTl and |Vl Figure [2] shows the measured value of T{p) as a function of the 
relative momentum squared. The amplitude T[p) is fitted with a naive polynomial function 
for m^ and p'^ 

A.oml + A2oml + A^oml + Aq^p^ + A^my , (38) 

where the indices of the parameter Aij denote the powers of ml and p"^, respectively. This fit 
form does not include p^ term, because our calculation is carried out with only two different 
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total momentum frames. In order to include a p'^ term in the fit form, we need several data 
with different relative momentum. Since the scattering amplitude at p = is nothing but 
the scattering length aom.,^, we need the ^30 term to coincide with the pion mass dependence 
of the scattering length in the previous section. The result of the fit parameters is presented 
in table IVIl and the fit lines for each m"^ are plotted in Fig. [21 

The ao/rriT^ estimated with the fit result, shown in Fig. [1] as denoted by "(5(p) analysis", 
is reasonably consistent with those of the CM analysis at each pion mass. We obtain 
o-o/'^TT = —1.99(12) GeV~^ at the physical pion mass, which agrees with the fit result 
given in the previous section (see table IIIII) . 

The measured scattering phase shift is plotted in Fig.Oas well as the result at the physical 
pion mass estimated from the T{p) fitting. The result is compared with the prediction of 



ChPT 581] estimated with experiment in the figure. 



3. Finite volume effect from final state interaction 

We estimate the derivative of the scattering phase shift from the fit results eq. fl38|) . and 
evaluate the derivatives for the functions 0(g) and (f)p{q) numerically to utilize the LL formula 
eq.([I]) and its extension eq.(I7]), respectively. 

For convenience, we define the conversion factors F and F^, 




9S dd) I .. 



p J { op uq ) 

which are the factors connecting the finite volume, decay amplitude in CM and Lab frames 
to the CM one in the infinite volume, respectively. The results of F and F^ are shown in 
tables [Til and |Vl To investigate the size of the interaction effect, we evaluate a ratio of the 
conversion factor to the one in the non-interacting case. We estimate the factors without 
the interaction F and F in each frame as, 

r = .IZAWh^, (42) 
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where the quantities with the overhne are the ones without two-pion interaction. In eq. (H2l) . 
Et^t^ and p are determined by 



t2 



K^ = A{ml+f) = {W'lf-P\ (43) 

with 



W^ = ^ml + P2 ± ^^. (44) 



In the CM frame case, p = 0, and then F = \J 4:{Et^t^LY as in Ref. 2J], while in the Lab 
frame case, 



f'' = J2WI {(wIy - {w'IY^ L\ (45) 



^p 



The ratios F/F and F^ / F are summarized in tables [TTl and |Vl The results show that 
the interaction effect is relatively large, e.g., 13 - 17%. The effect depends on the pion mass 
in the CM case, while there is no large dependence in the Lab case. In the CM case we find 
that the ratio decreases as the pion mass increases. This trend was also seen in the previous 



calculation 16| where the interaction effect was estimated by one-loop ChPT. It might be 
described by the fact that the scattering effect of the two-pion decreases in the smaller pion 
mass, because a^ oc m^. 

B. Decay amplitudes 

1. Off-shell decay amplitudes in finite volume 

To determine off-shell decay amplitudes in finite volume, we define ratios of correlation 
functions Ri{t) in the CM frame and Rf{t) in the Lab frame as 

T3P(f\ ^ /2 G'^ {t,tj,,tK)Z^^Zj^ 

'''^'' V3G^,(t,OG^(t,t,,)' ^^ 

where Gt,tt{G^^) and Gk{G^) are the 1 = 2 two-pion four-point function and the kaon two- 
point function in the CM (Lab) frame, respectively. These correlators are calculated with 
the wall sink operator as described in Sec. IlICi Gi and Gf are the A/ = 3/2 K ^ titi four- 
point functions defined in eqs. (iT3l) and (iT4l) . The index i denotes each operator i = 27,88 
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and m88 defined in eqs. (TT6|) . flT7|) . and flTSl) . zItt and Z\^ are the overlaps for tlie relevant 
operators witli eacli state. 

Wfien tlie correlation functions in eqs. fl46l) and WJ\ are dominated by eacli ground state, 
tfie ratios in principle will be a constant for those values of t, which corresponds to off-shell 
decay amplitude. We determine the desired amplitudes from the ratio in the flat region 
where the effective energy of each state has a plateau. The typical plateaus for the two-pion 
and kaon states are presented in Fig. |H 

The overlaps are determined by the wall sink correlators with the following fit forms, 

(^^^(tjt^) = Z^^ • ATTTri-E'TrTT) F ~ ^ttI); (48) 

GL^t.U) = (C)' ■ ^L{EL \t - U\), (49) 

GKit,tK) = Zl-AKimK,\t-tK\), (50) 

G^^{t,tK) = {zff-AK{E^,\t-tK\), (51) 

where E'^tt, E^^, ttik, and E^ are the measured energies from the point sink correlators, 
whose values are presented in tables. [TTl 13 IVIIl and IVIIIt respectively. The uik and E^ 
are in tx = 20 case. The kernels are defined by 

A^^{E,t) = e-^* + e-^(2^-*) + C, (52) 

A^„{E, t) = e-^* + e-^(2T-i) ^ ^ L-(i?.^— .)< + ^-{E^-m^){2T-t)'\ ^ (53) 

AK{E,t) = e-^* + e-^(2T-<)^ (54) 

where m^ and E^ in eq. flS51) are the measured single pion mass and energy, and C is one of 
the fit parameters. In the extraction of the overlaps, we use the fitting range of 6 < t < 25 
for the two-pion correlators, and < t < tx — 7 for the kaon correlators in both the frames. 

Figure [5] shows Ri{t) for all the operators with the different tx obtained in the CM 
calculation at the lightest 771^ and iris. All the ratios have clear signals, and are almost flat 
in the region oi t-^ <^t <^tK- Thus, the off-shell amplitudes are determined from averaged 
values in the fiat region of 6 < t < t/^ — 7. The values with the error are presented in each 
panel by the solid and dashed lines. The averaged values in different tx are consistent with 
each other. 

The ratios Rf{t) are shown in Fig. [6] at the same parameters as in the CM case. For all 
operators, the statistical error of the ratio is much larger than those in the CM case, while 
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the error decreases as tx decreases. In the case of tx = 25 it is hard to determine a flat 
region due to the huge error, while we can see a reasonable flat region for Ik = 16 and 20 
cases. We choose the same time slice range for the averaged values as in the CM case. The 
results of the off-shell amplitudes for tx = 16 and 20 reasonably agree with each other. 

According to the above discussion, the result in the t^ = 16 case has the smallest sta- 
tistical error in the Lab frame. The separation between the two-pion and the kaon sources, 
however, seems to be small in this case. Thus, we should worry about excited state contam- 
ination. In order to avoid this systematic error as much as possible, we choose the results 
with tx = 20 in the following analyses. We tabulate the result of the off-shell amphtudes 
for each 771^ and mg in tables llXt IX] and |Xl] for the CM calculation, and in tables IXIIt IXIIIt 
and IXI Vl for the Lab calculation in the case of Ik = 20. 

2. On-shell decay amplitudes in finite volume 

We determine the finite volume, on-shell, amplitudes in both the frames, \Mi\ and jM/'l, 
for each operator i = 27, 88 and m88, by interpolating the off-shell amplitudes with different 
strange quark masses at fixed m„. The on-shell kinematic point corresponds to rriK = E^,^ 
in the CM frame, and E^ = E^^ in the Lab frame. 

We plot the off-shell decay amplitude as a function of the kaon energy in Fig. [71 The 
figure shows that the amplitudes at rriu = 0.015 are described by a linear function of the kaon 
energy using three strange quark masses. On the other hand, the off-shell decay amplitudes 
for heavier m„ show large curvature with the kaon energy, so we employ a quadratic function 
of rriK or E^. The results of the on-shell amplitudes are presented by open symbols in the 
figures, and also tabulated in tables IS IS IXH IXHl IXTTTl and |XlV] denoted as "On-shell" 
for each operator and frame. 

3. On-shell decay amplitudes in infinite volume 

The on-shell decay amplitudes in infinite volume, \Ai\ for i = 27,88, and r7i88, are ob- 
tained for the LL formula eq. ([T]) and its extension eq. ([7]) by combining the measured on-shell 
amplitudes \Mi\ and \Mf\ in finite volume and, the conversion factors F and F^ in tables HTl 
and |Vl The results of the infinite volume, on-shell, decay amplitudes are summarized in 
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table IXVl 

We compare the infinite volume, on-shell decay amplitudes obtained from the different 
frame calculations. Figure M shows the decay amplitudes of the 27 operator with the CM 
and Lab calculations at the lightest m„ as a function of p^. We also plot previous results 



calculated with H-parity boundary condition 19|, |20| in the figure. The amplitude calculated 



in the Lab frame is consistent with the line interpolated between those of the CM and H- 
parity boundary condition calculations. This momentum dependence is consistent with a 
simple expectation that the Lab result is smoothly connected to the results obtained from 
the CM frame. In contrast to this result, the one obtained with tx = 16, also plotted in the 
figure, is below the linear fit. It may suggest that excited state contaminations are large in 
the tx = 16 case. 



C. Weak matrix elements 



The matching factors of the weak matrix elements were previously calculated in Ref. 



y 



using the regularization independent (RI) scheme and a non-perturbative method 63|, |64| 
at the scale n = 1.44 GeV. The renormalization factors are summarized below. 



/^ 0.8322(96) ^ 



Zij 



(55) 



0.8940(75) -0.0561(71) 

\ -0.0795(58) 0.964(17) / 

with the operator indices i,j = 27, 88,?7i88. In the following analysis we neglect the statis- 
tical errors of the renormalization constants, because those of the matrix element are larger 
than the ones of the diagonal parts. The renormalized weak matrix elements {Af^l = Zij\Aj\ 
for i,j = 27, 88 and r7i88 are summarized in table IXVll 

1. 21 operator 

The weak matrix element of the 27 operator is shown in Fig. [91 The horizontal axis is the 
pion mass squared. The matrix elements obtained from the CM and Lab calculations are 
clearly distinguished from each other. This means that the relative momentum dependence 
is important in the matrix element. The figure also shows that the matrix element decreases 
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with the pion mass. In order to investigate the dependences of the pion mass and relative 
momentum, we carry out a global fit of the result with a naive polynomial form, 

Boo + Bionil + Boip". (56) 

We omit a p^ term, because our result is obtained at only two different momentum in each 
pion mass, as discussed in Sec. IIVAI The value of x^/d.o.f. is 8.6 when we fit all the four 
pion mass data, so that we exclude the heaviest data to make the x^/d.o.f. acceptable. We 
obtain Boo = —0.0022(23) which is consistent with zero within the statistical error. This 
implies that the matrix element at p = vanishes in the chiral limit as shown in Fig. [9] as a 
dotted line. Thus, a fit form without the constant term Boo also gives a reasonable x^/d.o.f. 
as summarized in table IXVIII 

This tendency is consistent with the prediction of leading order (LO) ChPT y, |65|, 



^^rl = -«27^«-m2) (57) 



= -«27^(3m^ + V), (58) 

where a27 is a constant and / is the pion decay constant in the chiral limit. In the second 
line we use the on-shell condition, 

m\ = A{ml+p'), (59) 

because our matrix element is calculated for on-shell states. The normalization of eg. (1571) is 



based on Ref. ll|]. It should be noted that while the behavior of the matrix element near 
the chiral limit agrees with the LO ChPT, we omit logarithms predicted at NLO due to the 
large pion mass used in our calculation. 

We attempt to estimate a systematic error stemming from neglecting the NLO log terms 



in (quenched) ChPT formula [66(]. Both the formulae are explained in the appendix, and the 
results are summarized in table IXVIIII In the fits we assume that the formulae are valid in 
the lighter pion mass region, 771^^ < 0.5 GeV, and then we use only the data at the two lightest 
pion masses. In this restricted pion mass region, both the formulae work well. We estimate 
a dimensionless quantity 0:27//^ with / = 0.133 GeV [35| using the full (quenched) ChPT 
formula, and then find that the absolute value is 52(37)% smaller than the one estimated by 
Bio of the linear fit. This suggests that the determination of 0:27//^ has the large systematic 
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error. In order to remove the systematic error, we need to calculate more data at the region 
where ChPT is valid. While 027//^ varies largely with the fit form, the linear fit and the 
two ChPT fits give consistent results at the physical point (see tables IXVlII and [XVIIip . 



m^ = 140 MeV, p = 206 MeV, (60) 

which is determined by the on-shell condition with the physical kaon mass rriK = 498 MeV. 

2. 88 and m88 operators 

The weak matrix elements of the 88 and m88 operators are presented in Figs. [10] and 
im respectively. The matrix elements of the 88 operator obtained from the different frames 
are clearly separated and have appreciable slopes. The pion mass and relative momentum 
dependences are important in the matrix element as in the 27 operator case. In contrast 
to the 88 operator, the results of the m88 operator have neither strong pion mass nor 
relative momentum dependences. Both the matrix elements do not have clear curvatures 
with respect to the pion mass squared. Thus, we simply adopt the polynomial function 
in eq. flS^ for the global fit of the matrix elements. In this fit we use the data at all four 
pion masses. The fit results are tabulated in table IXVIIi In the m88 case the x^/d.o.f. is 
relatively larger than the other cases. This might be due to the poor determination of the 
covariance matrix of the matrix element. The matrix element at the physical point and the 
chiral limit is almost independent of fit form, because it does not have strong dependences 
for m'^ and p^. Therefore we choose to quote this fit as our result. 

We plot the fit results at the chiral limit with p = for the 88 and m88 operators in 
Figs. [To] and [Til respectively. A constant remains in both the matrix elements at the limits 
rTT-Tr = p = 0. While the trend is quite different from the 27 case, it is consistent with the 



prediction of LO ChPT |65l.l67. 



68 



69| 



Af\ = -ai^^, i = 88, m88, (61) 



where a^ is a constant. Again, although these trends agree with the behavior in the chiral 
limit of LO ChPT, we omit logarithms which enter NLO in ChPT due to large pion mass. 
A systematic error stemming from the fit without log terms is estimated using full ChPT 
formula f66|, (see Appendix for details). The quenched formula does not work in both the 
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matrix elements as presented in table IXVIIi Therefore, we will not discuss the fit result 
obtained with the quenched formula. We use the same assumption for the pion mass as in 
the 27 operator case, and then fit with only the two lightest pion mass data. It is found that 
a dimensionless quantity Oj//^ {i = 88, r?7,88) in the ChPT fit is reasonably consistent with 
that of the linear fit, determined by Bqq. The differences are less than 11%. This means that 
the low energy constants are less sensitive to the NLO logs than that in the 27 operator. As 
presented in tables IXVllI and IXVlllt the polynomial and ChPT fit results at the physical 
point, TTijr = 140 MeV and p = 206 MeV, also reasonably agree with each other. 

3. Comparison between direct and indirect methods 

We attempt to compare the results for the direct and indirect methods through the 
dimensionless parameters in LO ChPT, 027//^ and Oj//^ {i = 88,m88), which have been 
already discussed in the previous sections. The results are tabulated in table IXIXI as well 



as those obtained from the previous indirect calculation HI] . The constants of the previous 



work are estimated with / = 0.137 GeV [36| and the renormalization scale /i = 2.13 GeV. 
We find that the results for ass and ctmss are almost consistent in both the methods, while 
0:27 differs by a factor of 2.6. Our result of 0:27, however, contains a large systematic error 
of the chiral extrapolation, and varies by a factor of 1/2 as can be found in table IXIXi 
Other possible systematic errors in this comparison are different choice of parameters, such 
as lattice spacing and the renormalization scale, and different choice of gauge action. 

These parameters in LO ChPT obtained from the direct method should be consistent 
with those from the indirect method if no systematic error is included. This is because the 
final state interaction effect of the two pions in the matrix elements vanishes in both the 
limits m^ = p^ = 0. However, more strict consistency check of these parameters is beyond 
this work, because many systematic errors are included in the comparison discussed in the 
above. For this check, we need to calculate the matrix elements with both the direct and 
indirect methods on exactly same configurations at the lighter pion masses. 
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D. KeA2 and Im742 



We calculate IieA2 and I111A2 from the weak matrix elements l^^yl, l^^sl, and l^mssl- 
The definition of the decay amplitude Aj is given by 



{inn)i\Hw\K'') = Aje''' 
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^(2i(/x) + Tyi{fi)){Q,)i{fi) 



i=l 



(62) 
(63) 



where the index I denotes the isospin, Zi{fi) and yj(/i) are the Wilson coefficients at the scale 
H, and other parameters are presented in table IXX[ The relations between the weak matrix 
elements {Qi)2 and \Aj^\ for j = 27, 88,m88 are listed below, 



1^27 1 



\A 



RIi 



\A 



RI I 
mSSi 



3|(Ql)2h 

21(^7)2! 

2\{Q8)2\. 



S\{Q2)2\=2\{Q,)2\=2\{Q,o), 



(64) 
(65) 
(66) 



The other weak matrix elements vanish in the 1 = 2 case, i.e., ((53)2 = (<54)2 = (Q^)-?. 



20|. The 



(^6)2 = 0. The detail of the calculation of the Wilson coefficients is in Ref. 
coefficients are evaluated in the NDR scheme TO], which are converted to the RI scheme 
at the scale yU = 1.44 GeV [20|. The values of the Wilson coefficients are summarized in 
table IXXn 



1. ReA2 



The real part of A2 is evaluated through the following equation, 

10 



ReA. 



dp 

V2 



KJIK 



ud\ ]'' us\ 



^(^,(/x) + Re(r)y,(/x))|(Q,)2|(/i) 



1=1 



(67) 



Rey42 obtained in the CM and Lab frame calculations is tabulated in table IXXIII 

The upper panel of Fig. [12] shows that the main contribution of YieA2 comes from | {Qi)2\ 

and I ((52)21 in the CM calculation with the lightest pion mass (contributions come only from 

1^27!). This is consistent with previous results 10, lUl- The trend is not changed in the Lab 

calculation as shown in the lower panel of Fig. [121 

We will not discuss the chiral extrapolations of IieA2 in detail, because it is essentially 

same as those of | A^j \ in Sec. IIVC II In spite of this, we just summarize the result of the chiral 
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extrapolations in table IXXIIII At the physical point, we obtain Rey42 = 1.66(23) ( 



+48 N 
-03 J 



XlO" 



GeV, where the first and second errors are statistical and systematic, respectively. The 
central value and the first error are obtained from the simple polynomial fit eq. fl56l) with 
-Boo = 0. The systematic error is estimated by comparing the central value with the result 
from a fit with an added rn^p^ term. 



BiQinl + BqiP^ + Bumlp^, 



(68) 



and also (quenched) ChPT fit using the two lightest pion mass data as in Sec. IIV C 1[ 
The larger systematic error stems from the polynomial fit result with eq. ( I68ll as shown in 
table IXXIIH while the ChPT fits give consistent results with the one obtained from the 
simple polynomial fit, eq. fl56|) with Bqq = 0, as in \Afj\. 



In Fig. [13] we plot the resu 



t at the physical point as well as the previous results obtained 
with the indirect method 10, lUl, direct calculation with ChPT 16|, and the experiment. 
Our result reasonably agrees with the experimental value and also the previous results 
except the result in Ref. HI] (choice 2 with fi = 2.13 GeV). The difference of R.eA2 between 
our result and the RBC result is a factor of 1.4. This difference is smaller than what we 
observed in 027//"^ as described in Sec. IIVC II We may consider that ReA2 oc a27{3ml+4:p^) 
according to the LO ChPT relation eq. fl58l) . so that the reduction of the discrepancy between 
our result and the RBC result, found in 0:27//^ to that of R.eA2, seems inconsistent with 
the expectation derived from the LO relation. Indeed, our polynomial fit with Bqq = is 
not consistent with LO ChPT, as -Bio/3 > Bqi/4 can be read off in table IXXIIII although 
both the quantities give the same value of 0(27/ f^ if LO ChPT is valid. This inconsistency 
between Biq/3 and -Boi/4 causes the smaller difference between ours and the RBC result of 
ReA2 than that of a27/ f^- Again, this is due to the large systematic uncertainty during the 
determination of 0.211 f^- 



2. ImA2 

The imaginary part of A2 is determined through the definition of the decay amplitude, 
eq. (l63l) . as 



Imy42 = — =|Kd||Ks| 
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5^Im(r)i/,(/.)|(Q,)2|(/i) 



i=7 



(69) 



The results at each pion mass and frame are tabulated in table IXXIVI 
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Figure [m shows the contribution of each operator to Iniy42 in the CM and Lab calculations 
at the lightest pion mass. In both frames the largest contribution to Imy42 is \{Qs)2\, which 
is constructed only from |y4^gg|, and the second largest one is | ((59)21 with the opposite sign 
and comes from \A^j\. Figure [T^ shows that lmy42 has a significant slope for m^, while the 
leading contribution of ImA2, l^mssl does not have such a slope as can be seen in Fig. [TTl 
This slope is caused by the large pion mass dependence of |((59)2| and |(Qio)2|7 or \Afj\, as 
shown in Fig. [16] for the CM case. This tendency does not change in the Lab calculation. 

We estimate ITS1A2 at the physical point using the same polynomial fit form, eq. (!56l) . 
because the largest contribution of lm.A2 is given by l^mssl- The result with the phys- 
ical momentum and its chiral extrapolation are plotted in Fig. [T51 We obtain lm.A2 = 
— 1.181 (26) (^0^4) X 10^"*^^ GeV at the physical point. Again, the central value and statis- 
tical error are obtained from a fit to eq. fl56l) . while the systematic error is determined by 
comparing the central value with that from a fit form with an added m'^p'^ term. 

Boo + B.oml + Boip" + Bumlp\ (70) 

and ChPT formula for 88 and ?7i88 operators, explained in appendix, with the data at the 
two lightest pion masses as in Sec. IIVC2I In this case, the quenched ChPT formula does 
not give a reasonable x^/d.o.f. as in l^^gmssl- The larger systematic error comes from the 
fit with the ChPT formula. These fit results are presented in table IXXVI 

The result at the physical point is compared with the previous indirect calculation re- 
11| in Fig. [T71 The differences of our result from CP-PACS 72| and RBC (choice 
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2 with /i = 2.13 GeV) results are a factor of 1.9 and 0.93, respectively. However, from the 
comparison we cannot conclude that the final state interaction effect is appreciable or not in 
Im742, because the previous works employed different simulation parameters from our cal- 
culation, and also even the two results with the indirect method are inconsistent. Moreover 
our result has about 10% systematic error of the chiral extrapolation estimated in the above. 
More a detailed comparison between the direct and indirect methods is required for making 
a firm conclusion how much the final state interaction effect is in this quantity. 
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3. Systematic error from definition of momentum 

A systematic error for 'ReA2 and IniA2 arises from the violation of the naive dispersion 
relation on lattice. This error would decrease toward the continuum limit. To estimate the 
systematic error, we compare the results in the above Sees. IIVD II and 



22 



56| 



V D 21 with the ones 



given in the analysis with the lattice dispersion relation in Refs. 

cosh(Ef ) = cosh(m^) + 2 sm^{P/2). (71) 

In the analysis eqs.([2]) and ([9]) are replaced by 

2sin2(p/2) = cosh(K^/2) -cosh(m^), (72) 

cosh(E^^) = cosh(Efj - 2sin2(P/2), (73) 

respectively. Basically, these equations provide a different relative momentum p from the 
naive analysis. 

In the CM results, the change of the relative momentum is small, 0.6-1.8% in all the 
four masses, which increases as m.,; increases. The changes for Rey42 and Imy42 are less than 
1.5% in this case. In contrast to the CM case, p differs by 2.7-5.2% in the Lab calculation. 
The difference also increases as 171^^ increases. At the heaviest pion mass the analysis gives 
12(11)% smaller (larger) value than the naive dispersion analysis in 'ReA2{hivA2) . The main 
source of this difference is the conversion factor whose difference is 8.6% at the heaviest point. 
At the lightest pion mass the difference in Rey42(Imy42) is reduced to 5.8(5.2)%, where the 
difference of the conversion factor is 3.0%. This suggests that the systematic error decreases 
as the pion mass gets lighter at fixed relative momentum, or as the momentum decreases at 
fixed pion mass. 

In this analysis we determine IieA2 and Imy42 at the physical point with the same fit 
forms as used in the above sections. For Rey42 the fit form, eq. fl56l) with Bqq = 0, no longer 
gives a reasonable x^/d.o.f., because the Lab results at the heavier point are largely changed. 
On the other hand, eq. fl68|) gives a consistent result, Rey42 = 2.19(36) x 10"*^ GeV, at the 
physical point compared to the one given by using the same form as in Sec. IIVD 1[ For 
Im742 the consistent result is obtained with the fit form eq. flSBl) . while the result with eq. fITU]) 
is larger, lTaA2 = —1.137(37) x 10^^^ GeV, than the corresponding one in Sec. IIV D 21 

We estimate the systematic errors for IieA2 and lm.A2 at the physical point arising from 
this analysis by comparing the central values with the ones obtained in the above sections. 
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We combine them with the previous systematic errors, and finally quote 

ReA2 = 1.66(23) (+^^)(+f) x IQ-' GeV, (74) 

ImA, = -1.181(26)(+Jtl)(^o') x 10"^' GeV, (75) 

where the first error is statistic, and the second and third are systematic ones. 

V. CONCLUSIONS 

In this article we have presented our results of the A/ = 3/2 kaon weak matrix elements 
calculated with non-zero total momentum using the quenched approximation on a coarse 
lattice; a~^ = 1.31(4) GeV. The calculation is carried out with an extension of the Lellouch 
and Liischer formula, recently proposed by two groups, to obtain the infinite volume, on- 
shell, decay amplitude. It is very encouraging that we have obtained the on-shell weak matrix 
elements, taking into account final state interactions properly, with a reasonable statistical 
error. We have found that our result of R,eA2 at the physical point is reasonably consistent 
with the experimental value, and Im.A2 is comparable with previous results using the indirect 
method. While we have attempted to compare our results with those from the indirect 
method, many systematic errors are involved in the comparison, e.g., chiral extrapolations 
and different simulation parameters from previous calculations. A more comprehensive 
investigation is required. 

Another systematic uncertainty may come from the fact that we use the heavy strange 
quark masses, which we have used in order to evaluate the decay amplitude closely at the 
on-shell kinematical point in this study. The difference of the definition of the relative mo- 
mentum on lattice may cause an additional systematic error. To get rid of all the systematic 
errors related to a rigorous satisfaction of the on-shell condition, we need either the lighter 
pion mass or larger volume simulations to make the two-pion energy closer to the physical 
kaon mass, and also the simulations at finer lattice spacing are preferable. 

We also note that although chiral log behavior in the weak matrix elements are expected 
from ChPT, we have not seen such effects in our data. The problem may stem from the 
large pion mass and the coarse lattice spacing in our simulation because the prediction of 
ChPT is valid only in small pion mass region and may also be modified by non-zero lattice 
spacing effects. To confirm expected ChPT behavior, the simulations at the smaller pion 
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mass and finer lattice spacing are again required. 

Apparently the systematic error due to quenching is uncontrolled in this study, so that 
we should extend the present calculation in the dynamical lattice simulation as well. Besides 
the investigation of the unrevealed systematic errors, further investigation of the final state 
interaction remains an important future work. It is also important to calculate the A/ = 1/2 
kaon weak matrix elements with this method to evaluate the CP violation parameter e' je 
and the A/ = 1/2 selection rule. However, it is more difficult than the present A/ = 3/2 
case because the / = tttt scattering contribution enters in the A/ = 1/2 case. 
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APPENDIX A: NLO CHPT FORMULAE AT PHYSICAL KINEMATICS 



The NLO ChPT formula for the 27 and 88 operators at the physical kinematics in full 
and quenched case is obtained in Ref . 66| . Following the analysis in Ref . 2l|] , we rewrite /tt 
and Jk in the formula in terms of /. Thus, we use the following fit forms in ChPT analysis. 
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for the 88 (m88) operator, where Izf,a,b,c+d and Ja,b,c+d are presented in Ref. [66|, and It,^k 
in Ref. [TlJ. /j and Jj contain log terms and the scale fi. In our analysis we set /x = 1 GeV 
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for simplicity. We omit a p^ term in eq.OAlD as explained in Sec. II VC II The normalization 



of these equations are based on Ref. |ll|. 



The quenched formula is obtained by replacing hf.a.h.r-^d and Ja,b,c+d by /^j „ ;, c+d ^^'^ 
•^abc+d i^ ^^^- m], respectively, and Ik by //< in Ref. [7l|, and setting J^r = in the 



above equations. In the quenched analysis, we need two parameters a and niQ. We find 
that the result does not largely depend on the parameters in the ranges, < a < 0.5 and 
0.1 < mo < 0.866 GeV, so that we use a = 0.1 and mo = 0.866 GeV in our analysis. 
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[72] CP-PACS result is estimated from p(3/2) ^^ 24^ x 32 lattice of A^ = 325 MeV with a hnear 
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tx TTlu ''TT'S # of COnf. 

16, 20 0.015 0.12, 0.18, 0.24 371 

0.03, 0.04, 0.05 0.12, 0.18, 0.24, 0.28. 0.35, 0.44 252 

25 0.015, 0.03, 0.04, 0.05 0.12, 0.18, 0.24, 0.28. 0.35, 0.44 100 

TABLE I: Time slice of kaon operator t^, u,d quark mass m„, strange quark mass m^, and number 
of configuration. 

ruu 0.015 0.03 0.04 0.05 

m^[GeV] 0.35462(97) 0.4784(10) 0.54609(94) 0.60703(89) 

E^^[GeV] 0.7233(20) 0.9698(20) 1.1043(19) 1.2254(18) 

j9[GeV] 0.0710(19) 0.0792(16) 0.0816(16) 0.0831(16) 

6ip)[deg.] -2.55(19) -3.43(19) -3.72(19) -3.90(19) 

T{p) -0.227(11) -0.367(13) -0.439(14) -0.502(16) 

ao/m^[l/GeV2] -1.770(81) -1.580(54) -1.455(47) -1.348(42) 

q-id(j){q)/dq) 0.1215(83) 0.1601(81) 0.1724(82) 0.1801(83) 

p-{d6{p)/dp) -0.0473(32) -0.0611(29) -0.0657(30) -0.0685(31) 

F 44.41(48) 67.56(51) 81.50(54) 94.86(59) 

F/F 0.8708(89) 0.8456(59) 0.8363(52) 0.8306(49) 

TABLE II: Results of CM calculation. Pion mass m.j^, two-pion energy -Etttt, relative momentum p, 
scattering phase shift 6{p), scattering length over pion mass ao/rriT^, and derivatives q- {d(t)p{q)/dq) 
and p ■ {d6{p) / dp) are summarized. T{p) and F are scattering amplitude and conversion factor, 
which are defined in eqs.(|37p and (j39p . F/F is ratio of conversion factor to one in non- interacting 
case. Definition of F is given in eq. (j4ip . 



A[l/GeV2] B C[GeV2] xVd-o.f. mf^' 
-2.10(13) 2.78(69) -1.93(92) 0.11 -2.05(11) 

TABLE III: Fit result of scattering length aQ/m-,^ [1/GeV^] with a quadratic function. Results at 
physical pion mass are also listed. 
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/[GeV] /(/.) Q xVd.o.f. mf^^ 

0.1330(56) 0.85(14) 0.68(23) 0.22 -2.13(15) 

TABLE IV: Fit result of scattering length ao/m^ [l/GeV^] with NLO ChPT formula. Results at 
physical pion mass are also listed. 



niu 0.015 0.03 0.04 0.05 

m^[GeV] 0.35462(97) 0.4784(10) 0.54609(94) 0.60703(89) 

AE^^[GeY] 0.0250(36) 0.0213(20) 0.0196(15) 0.0183(12) 

£;^^[GeV] 1.0045(39) 1.2022(26) 1.3159(22) 1.4210(20) 

p[GeV] 0.2456(37) 0.2575(24) 0.2618(19) 0.2649(16) 

7 1.1643(16) 1.10641(54) 1.08646(33) 1.07277(23) 

5{p)[deg.] -11.6(1.6) -12.4(1.1) -12.59(92) -12.78(79) 

T{p) -0.362(48) -0.462(40) -0.517(36) -0.567(33) 

q-{d(j)p{q)/dq) 1.878(32) 2.137(20) 2.232(16) 2.300(13) 

p-{d5{p)/dp) -0.261(57) -0.266(42) -0.261(35) -0.253(30) 

F^ 48.9(1.2) 65.8(1.1) 76.1(1.0) 86.0(1.0) 

F^/F^ 0.869(20) 0.865(14) 0.865(11) 0.866(10) 

TABLE V: Results of Lab calculation. Pion mass m.,^, energy shift AE^^, two-pion energy E^^, 
relative momentum p, boost factor 7, scattering phase shift 6{p), and derivatives q ■ id(pp{q)/dq) 
and p • {d5{p) / dp) are summarized. T{p) and F^ are scattering amplitude and conversion factor, 

which are defined in eqs. (|37p and (I40p . F /F is ratio of conversion factor to one in non- interacting 

P 

case. Definition of F is given in eq.(|42p. 



Aio[l/GeV2] 


yl2o[l/GeV4] 


A3o[l/GeV6] 


Aoi[l/GeV2] 


.4ii[l/GeV4] 


xVd.o.f. 


-2.04(13) 


2.59(71) 


-1.78(93) 


-2.42(81) 


3.1(1.8) 
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TABLE VI: Fit results of scattering amplitude T{p) with polynomial function eq.(p8 
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niK [GeV] rris 

niu 0.12 0.18 0.24 0.28 0.35 0.44 

0.015 0.7073(12) 0.8507(13) 0.9753(15) _ _ _ 

0.03 0.7428(12) 0.8814(12) 1.0029(13) 1.0763(14) 1.1914(16) 1.3163(18) 

0.04 0.7666(11) 0.9023(11) 1.0221(12) 1.0946(13) 1.2086(14) 1.3324(15) 

0.05 0.7901(10) 0.9232(11) 1.0413(11) 1.1129(12) 1.2258(12) 1.3487(13) 

TABLE VII: Kaon mass for each light and strange quark masses in tx = 20 case. 



E'^lGeV] 



m. 



rriu 0.12 0.18 0.24 0.28 0.35 0.44 

0.015 0.8673(25) 0.9847(24) 1.0914(25) _ _ _ 

0.03 0.8948(21) 1.0111(19) 1.1168(19) 1.1818(20) 1.2856(21) 1.4002(22) 

0.04 0.9146(18) 1.0293(17) 1.1339(16) 1.1984(17) 1.3014(17) 1.4152(18) 

0.05 0.9343(16) 1.0476(15) 1.1511(14) 1.2150(14) 1.3172(15) 1.4303(16) 

TABLE VIII: Kaon energy with momentum P = 2it/L for all light and strange quark masses in 
tx = 20 case. 



|M27|[10-3] ms 

TTiu 0.12 0.18 0.24 0.28 0.35 0.44 On-shell 

0.015 1.026(16) 1.082(18) 1.118(21) — — — 1.033(16) 

0.03 1.200(17) 1.232(18) 1.246(20) 1.245(21) 1.227(22) 1.174(24) 1.246(19) 

0.04 1.328(17) 1.349(18) 1.353(20) 1.347(20) 1.320(21) 1.255(22) 1.343(20) 

0.05 1.453(18) 1.465(19) 1.462(20) 1.451(21) 1.416(21) 1.340(22) 1.410(21) 

TABLE IX: Decay amplitude of 27 operator for all light and strange quark masses in CM frame. 
Result is in lattice unit. Results of on-shell amplitude are also included. 
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|M88|[10-3] m, 

m„ 0.12 0.18 0.24 0.28 0.35 0.44 On-shell 

0.015 7.00(12) 6.56(11) 6.22(12) _ _ _ 6.95(11) 

0.03 4.709(76) 4.391(75) 4.133(74) 3.972(74) 3.693(74) 3.320(72) 4.216(74) 

0.04 3.996(60) 3.721(58) 3.494(57) 3.353(56) 3.108(56) 2.786(54) 3.326(56) 

0.05 3.511(51) 3.268(48) 3.065(47) 2.937(46) 2.718(45) 2.430(43) 2.705(45) 

TABLE X: Decay amplitude of 88 operator for all light and strange quark masses in CM frame. 
Result is in lattice unit. Results of on-shell amplitude are also included. 



|M^88|[10-2] rus 

Tiiu 0.12 0.18 0.24 0.28 0.35 0.44 On-sheh 

0.015 2.548(40) 2.417(39) 2.314(40) _ _ _ 2.532(39) 

0.03 1.837(28) 1.740(27) 1.660(27) 1.608(27) 1.513(27) 1.379(26) 1.687(27) 

0.04 1.632(23) 1.546(22) 1.473(21) 1.425(21) 1.339(21) 1.219(20) 1.415(21) 

0.05 1.499(20) 1.421(19) 1.354(18) 1.309(18) 1.229(18) 1.118(17) 1.224(17) 

TABLE XI: Decay amplitude of r7T,88 operator for all light and strange quark masses in CM frame. 
Result is in lattice unit. Results of on-shell amplitude are also included. 



|Mf7|[10-3] nis 

niu 0.12 0.18 0.24 0.28 0.35 0.44 On-shell 

0.015 1.31(11) 1.37(11) 1.41(12) — — — 1.38(11) 

0.03 1.500(70) 1.545(70) 1.562(72) 1.559(74) 1.531(76) 1.457(77) 1.554(74) 

0.04 1.630(56) 1.679(56) 1.696(57) 1.693(58) 1.663(59) 1.583(59) 1.651(59) 

0.05 1.774(49) 1.824(49) 1.841(49) 1.837(50) 1.802(51) 1.715(51) 1.726(51) 

TABLE XII: Decay amplitude of 27 operator for all light and strange quark masses in Lab frame. 
Result is in lattice unit. Results of on-shell amplitude are also included. 



35 



\M[s\[W-^] 

ruu 0.12 0.18 0.24 0.28 0.35 0.44 On-shell 

0.015 6.31(56) 6.19(48) 6.05(45) _ _ _ 6.16(48) 

0.03 4.41(24) 4.26(22) 4.09(21) 3.96(20) 3.72(20) 3.38(19) 3.92(20) 

0.04 3.85(16) 3.70(15) 3.54(14) 3.43(14) 3.21(13) 2.90(13) 3.16(13) 

0.05 3.47(12) 3.33(11) 3.18(10) 3.07(10) 2.870(97) 2.586(92) 2.614(93) 

TABLE XIII: Decay amplitude of 88 operator for all light and strange quark masses in Lab frame. 
Result is in lattice unit. Results of on-shell amplitude are also included. 



I^^88l[10-^] 

m„ 0.12 0.18 0.24 0.28 0.35 0.44 On-shell 

0.015 2.44(19) 2.42(16) 2.38(15) _ _ _ 2.41(16) 

0.03 1.843(84) 1.808(75) 1.757(70) 1.717(68) 1.632(65) 1.503(62) 1.701(67) 

0.04 1.688(58) 1.652(51) 1.603(48) 1.564(47) 1.484(45) 1.363(43) 1.467(44) 

0.05 1.595(45) 1.559(40) 1.511(37) 1.473(36) 1.396(35) 1.280(33) 1.292(33) 

TABLE XIV: Decay amplitude of m88 operator for all light and strange quark masses in Lab 
frame. Result is in lattice unit. Results of on-shell amplitude are also included. 



CM Lab 

mu |^27|[GeV3] |^88|[GeV3] \AmS8\[GeY^] \A27\[Ge\^] \As8\[Ge\^] \Am88\[GeY^] 

0.015 0.1031(20) 0.694(14) 2.528(48) 0.151(12) 0.676(53) 2.65(18) 

0.03 0.1892(33) 0.640(13) 2.562(47) 0.230(12) 0.579(30) 2.51(10) 

0.04 0.2461(41) 0.609(11) 2.593(44) 0.283(11) 0.541(23) 2.510(83) 

0.05 0.3007(50) 0.577(10) 2.609(41) 0.334(11) 0.506(19) 2.499(71) 

TABLE XV: Decay amplitudes in infinite volume for 27, 88, m88 operators in CM and Lab calcu- 
lations. 
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CM Lab 

m. \A^I\[GeY^] \Afi\[GeY'] \A%s\[GeY'] l-^WlJGeV^] |^fi|[GeV=^] lA^^sslJGeV^] 

0.015 0.0858(17) 0.4782(96) 2.382(46) 0.126(10) 0.456(38) 2.50(17) 

0.03 0.1575(28) 0.4288(88) 2.419(44) 0.1913(97) 0.377(22) 2.378(98) 

0.04 0.2048(34) 0.3992(78) 2.451(41) 0.2351(90) 0.343(17) 2.376(78) 

0.05 0.2502(41) 0.3693(70) 2.470(39) 0.2778(89) 0.312(14) 2.369(67) 

TABLE XVI: Weak matrix elements in regularization independent scheme for 27, 88, rn-88 operators 
in CM and Lab calculations. 



opr. Boo[GeY^] -Bio[GeV] 5oi[GeV] xVd.o.f. phys.[GeV3] 

27 — 0.6687(95) 0.48(11) 1.26 0.0336(47) 

88 0.534(11) -0.429(22) -0.86(12) 0.22 0.490(11) 

m88 2.391(49) 0.14(10) -1.89(63) 2.48 2.313(55) 

TABLE XVIL Fit results for weak matrix elements of each operator in regularization independent 
scheme, l^f' | for i = 27, 88, ?ti88. Fit function is defined in eq. (|56p . Results at physical point, 
m^ = 140 MeV and p = 206 MeV, are also tabulated. 



opr. 


«27//" /32o[l/GeV2] /3n[l/GeV2] xVd-o.f. 


phys.[GeV3] 


27(Full) 
27(Quen.) 


-0.0391(15) -20.4(6.3) 8.5(1.5) 0.02 
-0.0509(20) -11.7(4.8) 5.2(1.2) 0.08 


0.0331(13) 
0.0374(15) 


opr. 


ai/f 7io[l/GeV2] 7oi[l/GeV2] xVd-o.f. 


phys.[GeV3] 



88(Full) 


-5.17(19) 


-0.51(69) 


0.04(16) 


1.44 


0.432(20) 


88(Quen.) 


-8.20(31) 


0.28(44) 


-0.30(10) 


5.39 


0.554(24) 


m88(Fun) 


-22.37(90) 


-0.78(70) 


0.38(16) 


2.33 


2.039(94) 


m88(Quen.) 


-35.5(1.5) 


0.15(44) 


-0.09(10) 


7.24 


2.56(11) 



TABLE XVIIL Fit results with full and quenched NLO ChPT formulae for weak matrix elements 
of each operator in regularization independent scheme, \Af-^\ for i = 27, 88,m88. Fit functions 
are defined in appendix. Results at physical point, itIt^ = 140 MeV and p = 206 MeV, are also 
tabulated. 
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g-MGeV] /[GeV] a^i/f ass/f amss/f 

direct(Lin.) 1.31 0.133 -0.0809(11) -5.52(11) -24.68(51) 

direct(ChPT) 1.31 0.133 -0.0391(15) -5.17(19) -22.37(90) 

indirect 1.92 0.137 -0.0306(13) -5.89(30) -21.4(1.2) 

TABLE XIX: Comparison of dimensionless constant of LO ChPT. This work and previous work 11 1 
are denoted by "direct" and "indirect", respectively. 



Gf 1.166x10-5 GeV-^ 
\Vus\ 0.2237 

\Vud\ 0.9747 

l^i^l 0.0410 

Vtd 0.00708-0.00297i 
r 0.00133-0.000559i 



TABLE XX: Parameters for weak matrix element 



m 
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z-iifJ-) 



ViifJ-) 



1 


-0.3522 


0.0 


2 


1.17721 


0.0 


3 


0.00446831 


0.0241094 


4 


-0.0140925 


-0.0503954 


5 


0.00506909 


0.00563178 


6 


-0.015967 


-0.0928098 


7 


0.0000502692 


-0.000186283 


8 


-0.0000134347 


0.00118057 


9 


0.0000428969 


-0.0114749 


10 


0.0000117198 


0.0037748 



TABLE XXI: Results of Wilson coefficients 



2C|. 



m„ 



0.015 



0.03 



0.04 



0.05 



CM 4.241(81) 7.78(14) 10.12(17) 12.37(20) 
Lab 6.22(50) 9.45(48) 11.62(45) 13.74(44) 

TABLE XXII: Results of Re^2[10"^GeV] in CM and Lab calculations. 



Poly. Sio[10-VGeV] 


5oi[10-7GeV] 


Bu[W''' /GeV^] 


xVd.o.f. phys.[10-«GeV] 


3.306(47) 
3.304(47) 


2.39(56) 
3.60(99) 


-4.0(2.7) 


1.26 1.66(23) 
0.90 2.14(40) 


ChPT 12V3a^^^Jf[W-yGeY] 


/?2o[l/GeV2] 


/3n[l/GeV2] 


xVd.o.f. phys.[10-SGeV] 


Full -5.33(21) 
Quen. -6.93(27) 


-20.4(6.3) 
-11.7(4.8) 


8.5(1.5) 
5.2(1.2) 


0.02 1.636(65) 
0.08 1.847(72) 



TABLE XXIII: Fit results with polynomial and NLO full (quenched) ChPT formula of ReA2, and 
result at physical point, rriTj = 140 MeV and p = 206 MeV. 
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niu 0.015 0.03 0.04 0.05 

CM -1.036(22) -0.786(19) -0.624(17) -0.462(16) 
Lab -0.952(90) -0.635(53) -0.468(42) -0.301(35) 

TABLE XXIV: Results of lm.A2[l^-^^G(N] in CM and Lab calculations. 



Poly. 5oo[10-i2GeV] ^lollO-^VGeV] 


5oi[10-iVGeV] 


5n[10-^VGeV3] 


xVd.o.f. phys.[10-i2GeV] 


-1.333(25) 2.307(56) 
-1.332(25) 2.302(57) 


2.54(36) 
2.2(1.0) 


0.7(1.9) 


0.15 -1.181(26) 
0.15 -1.194(43) 


ChPT 24V3a,^^^//3[10-i2GeV] 


7io[l/GeV2] 


7oi[l/GeV2] 


xVd.o.f. phys.[10-i2QgY] 


1.250(42) 


-1.24(66) 


-0.02(16) 


1.24 -1.040(46) 



TABLE XXV: Fit results with polynomial and NLO ChPT formula of ImA2 , and result at physical 
point, m^ = 140 MeV and p = 206 MeV. 
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FIG. 1: Scattering length over pion mass and its chiral extrapolation. Star and cross symbols are 
prediction of ChPT and experimental result, respectively. Open circles are obtained from analysis 
of scattering phase shift, and are slightly shifted to minus direction in x-axis. 
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CM(Lab) calculation. 
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error, respectively. 
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FIG. 7: Interpolations of ofF-shell decay amplitude |Mj| and \M[\ for 27, 88,m88 operators to 
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FIG. 8: Comparison of decay amplitudes obtained from different calculations in infinite volume 
of 27 operator at ra^ = 0.015. 
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FIG. 9: Weak matrix element of 27 operator. Dotted line is obtained from global fit with p = 0. 
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FIG. 10: Weak matrix element of 88 operator. Dotted line is obtained from global fit with at 
p = 0. Dashed line and triangle symbol denote global fit result with physical momentum p = 206 
MeV. Square symbol is slightly shifted to minus direction in x-axis. 
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FIG. 11: Weak matrix element of m88 operator. Dotted line is obtained from global fit with 
p = 0. Dashed line and triangle symbol denote global fit result with physical momentum p = 206 
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FIG. 12: Absolute values of each weak operator contribution for Re^2 at lightest quark mass in 
CM(upper panel) and Lab(lower panel) calculations. Open and stripe bars denote positive and 
negative contributions, respectively. Signs of operator 1 and 8 contributions are changed. 
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FIG. 14: Absolute values of each weak operator Qi contribution for Im^2 at lightest quark mass 
in CM(upper panel) and Lab(lower panel) calculations. Open and stripe bars denote positive and 
negative contributions, respectively. Signs of operator 8 and 10 contributions are changed. 
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FIG. 15: 1111^2 obtained from CM and Lab calculations and its global fit. Triangle symbol denotes 
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